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Thermodynamics of topological nonlinear charged Lifshitz black holes 
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In this paper, we construct a new class of analytic topological Lifshitz black holes with constant 
curvature horizon in the presence of power-law Maxwell field in four and higher dimensions. We find 
that in order to obtain these exact Lifshitz solutions, we need a dilaton and at least three electro¬ 
magnetic fields. Interestingly enough, we find that the reality of the charge of the electromagnetic 
field which is needed for having solutions with curved horizon rules out black holes with hyperbolic 
horizon. Next, we study the thermodynamics of these nonlinear charged Lifshitz black holes with 
spherical and flat horizons by calculating all the conserved and thermodynamic quantities of the 
solutions. Furthermore, we obtain a generalized Smarr formula and show that the first law of ther¬ 
modynamics is satisfied. We also perform a stability analysis in both canonical and grand-canonical 
ensemble. We find that the solutions are thermally stable in a proper ranges of the metric parame¬ 
ters. Finally, we comment on the dynamical stability of the obtained solutions under perturbations 
in four dimensions. 
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I. INTRODUCTION 


The idea of correspondence between gravity in an anti- 
de Sitter spacetime and the conformal field theory living 
on its boundary (AdS/CFT correspondence), was first 
suggested by Maldacena in 1998 [1]. According to Mal- 
dacena’s arguments the effects of the string or M-theory 
in a d-dimensional AdS„+i x spacetime can be ap¬ 

peared in the form of a field theory on an n-dimensional 
r-constant brane which is the boundary of AdS„+i space- 
time. He expressed that the finite temperature configu¬ 
rations in the field theory which lies on the n-dimensional 
brane and is decoupled from the bulk correspond to black 
hole configurations in AdS„+i spacetime [1]. This pro¬ 
posal makes us able to study the non-perturbative as¬ 
pects of the field theories. This idea attracted a lot of in¬ 
terests rapidly and has been studied from different points 
of view [2] . Since the line element of the AdS spacetime, 

2 /2 n—1 

ds^ = —^dt^ + X] (^) 

is invariant under an isotropic conformal transformation 

t —y At, Xi —y \xi, T —y A T, (2) 

so the application of AdS/CFT is restricted to systems 
respected isotropic scale invariance. But, quantum criti¬ 
cal systems show scaling symmetry 


t —^ Xi —i Xxi, (3) 


‘Electronic address: mkzangeneh@shirazu.ac.ir 
^Electronic address: asheykhi®shirazu.ac.ir 
^Electronic address: mhd@shirazu.ac.ir 


where z > 1 is a dynamical critical exponent which shows 
the degrees of anisotropy between space and time. There¬ 
fore, we can generalize AdS spacetime to other space- 
times with anisotropic scale invariance 


ds^ = 

I2z 




+ r^Y.'^xl 


( 4 ) 


i=l 


This spacetime was first introduced in [3] and is called 
Lifshitz spacetime. Black hole configurations in the Lif¬ 
shitz spacetime are dual to the field theories enjoying 
anisotropic scale invariance on the boundary (brane). 

It is well-known that Lifshitz spacetime is not a vac¬ 
uum solution of Einstein gravity and it needs some mat¬ 
ter sources. In many works this matter source was con¬ 
sidered to be a massive gauge field [4-7]. Some of these 
solutions were obtained for fixed z [4-6]. For instance, 
the authors of [4, 5] focused on z = 2 case. The main dis¬ 
advantage of such models containing massive gauge fields 
is that it is usually impossible to obtain an analytic solu¬ 
tion for an arbitrary z, although some efforts have been 
made to do that [8]. In Ref. [9], it was shown that by 
considering a dilaton field, instead of a massive gauge 
field, we are able to derive exact solutions which helps 
us to get more insight. In addition, string theory in its 
low energy limit recovers Einstein gravity with a scalar 
dilaton field nonminimally coupled to gravity and other 
fields such as gauge fields [10]. Thermal behavior of un¬ 
charged Lifshitz black branes in Einstein-dilaton gravity 
has been investigated in [11]. Thermodynamics of lin¬ 
early charged Lifshitz black hole/branes in the context 
of Einstein-dilaton gravity has been studied in [9] . 

Matter sources possessing conformal invariance have 
been always of much interests [12]. The first black hole 
solution with a conformally invariance matter source is 
the Reissner-Nordstrbm solution in four dimensions. In 
higher dimensions, Maxwell source is no longer confor¬ 
mally invariant. However, as in the case of massless 




2 


Klein-Gordon Lagrangian, a particular power of Maxwell 
Lagrangian can respect conformal invariance in arbitrary 
dimensions [13]. In fact, in (n + l)-dimensional spacetime 
the action 


where 


S'pM = [ (5) 

J M 


^PM 




(ti+1)/4 


is conformally invariant, i.e., it is invariant under the 
conformal transformation —>■ Afj_ —>■ and 

$ —Black hole solutions with power-law Maxwell 
electrodynamic fields have been studied in [13-17]. 

Here, we shall consider a class of asymptotically Lif- 
shitz black hole solutions of Einstein-dilaton gravity with 
power-law Maxwell field. Although electrodynamic La¬ 
grangian is conformally invariant in (n-l-1)/4 dimensions, 
this does not implies that the power should be fixed for 
general cases. We first consider a four dimensional toy 
model and fully describe the procedure of solving held 
equations and gaining desired solutions. Then, we extend 
our solutions to the higher-dimensional spacetimes. From 
the prediction of string theory, we know that the space- 
time may have more than four dimensions. Therefore we 
get enough motivation for studying black hole solutions 
of Einstein gravity in all higher dimensions. Although 
it was thought for a while that the extra spacial dimen¬ 
sions are of the order of the Planck scale, recent theories 
suggest that if we live on a three dimensional brane em¬ 
bedded in a higher dimensional bulk, it is quite possible 
to have the extra dimensions relatively large and still 
unobservable [18, 19]. In this scenario, all gravitational 
objects including black holes are higher-dimensional. 

This paper is structured as follows. In the next sec¬ 
tion, we introduce the four dimensional action in the pres¬ 
ence of a general nonlinear electrodynamic and derive the 
equations of motion by varying the action. We also give 
some general remarks about nonlinearly charged topo¬ 
logical Lifshitz black holes. Furthermore, focusing on 
power-law nonlinear electrodynamic source, we compre¬ 
hensively describe the procedure of obtaining asymptotic 
Lifshitz topological black hole exact solutions. In section 
III, we generalize our study to higher dimensions. In 
section IV, we calculate the mass of black hole solutions 
by the modified Brown-York method and study thermo¬ 
dynamics of the obtained solutions. We also check the 
validity of the first law of thermodynamics by finding the 
Smarr formula for the mass. In section V, we investigate 
thermal stability of the solutions in both canonical and 
grand-canonical ensembles. The dynamical stability of 
4-dimensional AdS black holes is examined in VI. In the 
last section we give a summary and closing remarks on 
our work. 


II. ACTION AND ASYMPTOTIC LIFSHITZ 
SOLUTIONS 


The action of Einstein-dilaton gravity in the presence 
of a nonlinear electromagnetic and two linear Maxwell 
fields can be written as: 

-2A4 + /:(F,$)-^e-2^^^iL.l, (6) 

i=2 J 


where TZ is the Ricci scalar on manifold A4, 4> is the 
dilaton field, A 4 and A^’s are constant parameters. In 
Eq. ( 6 ) F and iJ^’s are the Maxwell invariants of elec¬ 
tromagnetic fields = d[f,A^] and 
where and {Bi)^s are the electromagnetic potentials. 
C{F, <I>) stands for nonlinear electromagnetic Lagrangian. 
Varying the action ( 6 ) with respect to the metric 
the dilaton field $, electromagnetic potentials and 
(Hi)^’s, lead to the following field equations 


pv = [ 2 A 4 - 1 - 2CfF — C{F, $) 


3 

3 




A 

i^2 

i^2 


2CfF^xF, ^ + 2 a^$a.$, 

(7) 

4 

4=9 

(8) 

{CFFf^n = 0, 

(9) 

= 0. 

(10) 


where we use the convention Xy = dXfdY. Our pur¬ 
pose here is to find asymptotic Lifshitz topological black 
hole solutions of the field equations (7)-(10). The line el¬ 
ements of such a metric in four dimensions can be written 
as [4, 9] 


= + ( 11 ) 

where 

( d0^ + sin^ {0)dip'^ k = 1 

dnl = } d0'^ + dcj)'^ k = 0 , ( 12 ) 

(_ d0^ + sinh'^{0)d(p k = —1 

represents a 2 -dimensional hypersurface with constant 
curvature 2k and volume lu. For k = 1, the topology 
of the event horizon is a two-sphere S^, and the space- 
time has the topology E? x . For fc = 0, the topology 
of the event horizon is that of a torus and the spacetime 
has the topology x T^. For /c = — 1, the surface S 






3 


is a 2-dimensional hypersurface with constant neg¬ 
ative curvature. In this case the topology of spacetime 
is X . It is worth mentioning that the asymptotic 
behavior of metric (11) should be 

kP 

/W = l + ^- (13) 

Note that this asymptotic behavior is the exact Lifshitz 
solution which is obtained in Ref. [4] for z = 2 and also 
it reduces to AdS solution for z = 1 as one expects. 

In this paper, we intend to consider the power-law 
Maxwell nonlinear source, but here we pause to present 
some general results coming from general nonlinear elec¬ 
tromagnetic sources. Using (11), (10) can immediately 
be integrated as 

(i/Ort = (14) 

Substituting the line element (11) and (14) to (8), one 
obtains 


r.3 + -2 




Al 


2 z 3 


p3 —2 / ^ 
2 = 2 


Kq, e 




(15) 


where the prime denotes derivative with respect to r. 
Let’s consider the case where A^’s are positive (later, 
we see that this is true) and £$ < 0 (in most of non¬ 
linear electromagnetic cases this is a true consideration 


[20]). In this case the RHS of (15) is always positive 
and therefore is an increasing function. On the 

other hand, we can prove that our black hole has just 
one horizon in this case. Consider that there are two 
inner and outer horizons (denoted by r_ and r_|_ respec¬ 
tively). In one side is an increasing function, 

while on the other side f _4r3+^$7 |^_= 0. 

This contradiction means that our assumption of having 
inner and outer horizon is not correct. Thus, if a black 
hole exists, it has just one horizon and its causal struc¬ 
ture is Schwarzschild-like. Since an increasing function 
has just one zero, if there is any, the zero of oc¬ 

curs at the horizon when we have black hole solutions. 
Consequently, < 0 inside the black hole horizon 

and > 0 outside of it and therefore since / < 0 

and / > 0, inside and outside the horizon of black hole 
respectively , we always have $' > 0. 

After presenting the above general results about a 
general nonlinear electromagnetic sources, we continue 
our study with power-law Maxwell source with the La- 
grangian £(U, $) = (—In this case, one can 
integrate the electromagnetic field equation (9) as 


Frt = (16) 

Substituting (11), (14) and (16) into Eqs. (7) and (8), 
one receives 


^/' +3/-br24>'2/ k (2p- ^ g2A,3.^2 

--fA4- — -f +2^7^ 


p 

rf + f + 2 fz-r^<^'^f 
P 


'2^i—p^z\\. — z)pp^p/\zp—i) X ^ j2 — 2zpA 

2=2 


2l-p;2(l-2)p^4p/(2p-l) 

k (2 p-l)g2 V2^iP$)/(2p-i) 


= 0, 


2 ^-PpP-PPr^P/PP-P 


3 9 \ • <f» 2 

y ‘ =0 

^2-22^4 


2=2 


(17) 

(18) 


r2/" + 3r(l+z)/' + 2r2$'2/ + 2(7+ Z+l)/ ^ g 2 pg( 2 Aip$)/( 2 p-l) 3 ^ 2 A .$^2 

-+A4- „, -2^ ,2_2, . =0^ (19) 


2 P 


2^-Pp (l-z)P7-4p/(2p-l) 

r2/$"+(r2f+ (3 + z)r/)$' X^pql^ 

P 


2=2 


__ ^ ^ 

Ol—P/S (1 —2)p^4p/(2p—1) j2—2zpA 

2=2 


( 20 ) 


r 


Subtracting (18) from (17), we arrive at 
(1 - z) -I- r2$'2 

with the solution 


function as 


where 6 is a constant and z > 1. With (22) in hand, 
by solving field equations (18)-(20), we can find metric 


= 0, 

(21) 

kP A 4 P m 

zr"^ z + 2 r^+2 

())■ 


2 P-ig+pp(2-i)+2 ^ 2 p — \Y ^(“2\/-i+*>‘ip)/(2p-i) 

(22) 

(2V-1 + Z\ip + Z{2p - 1) - 2) 7P(2-Ai V^)/(2p-l) 


-E 


qt 


,2^2 z j. —4+2 y/ —1+5;Ai 


^ (-2 + Z + 2V-1 + ZA) 62V^T+^aC 


(23) 
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The field equations are fully satisfied provided: 


and therefore the gauge potential is 


Ai 

A 2 

A 3 


-A 4 (z - 1 ) 54 
(z + l)Z2(^-i)’ 
kb"^ (z - 1 ) 
;2(z-1)2 ’ 

(1 - 2p) Vz-1 
P 

2 

\/z — 1 ’ 

1 

y/ Z — I 


q^fo 2 (z-l) 

r4rr4 


(29) 


As we will show later in section (IV), the dependence of 
total mass of the black hole on At is linear and there¬ 
fore in order to have a finite mass, r 4 should be posi¬ 
tive. Thus, the metric function (27) has the appropriate 
asymptotic given in Eq. (13). 


(24) 


A. PROPERTIES OE THE SOLUTIONS 


Using (24), the metric function /(r) may be written as 
kf 2 AaP 


fir) = 


m 


2z2 ( 2 : -I- 1) (2 -I- 2) r^+2 

2P-I^2p^2p(z-1)-H2 (2^ _ 1) 52 (z- 1 ) 
P4pr4+z-t-2 


(25) 


In order to study the properties of the solutions, we 
first mention the constraint on z and p which comes from 
the positivity of r 4 : 


z — 2 4 - 


2 

(2p-l) 


> 0 . 


(30) 


where r 4 = z — 2 -|- 2/{2p — 1). Looking at in (24), 
we find that the hyperbolic case k = — \ causes an imag¬ 
inary charge except for z = 1. Therefore, we continue 
the paper just for the cases k = 0 and k = 1. It is re¬ 
markable to note that our solutions include topological 
asymptotic AdS black holes for z = 1. At first glance the 
constants (24) are infinite for z = 1, but one should note 
that although A^s diverge at z = 1 , Aid) is finite for this 
case and = 0 and therefore action ( 6 ) reduces to the 
Einstein action in the presence of cosmological constant 
and nonlinear electromagnetic field for z = 1. This kind 
of solutions in the absence of cosmological constant has 
been introduced in [14]. Also, this kind of solutions in 
the presence of cosmological constant and in the context 
of Lovelock gravity has been studied in [16]. 

In order to have asymptotic Lifshitz solutions, /(r) 
should reduce to the metric (13). Thus, we have another 
restriction on A 4 as 


The inequality (30) imposes the following restrictions on 
p and z, 

for p< 1/2, z-1 > (3-2p)/(l-2p), 

for 1/2 < p < 3/2, all z(> 1) values are allowed, (31) 

for p> 3/2, z — 1 > (2p — 3)/(2p — 1). 

Eor the allowed ranges of p and z, the term including q in 
/(r) is dominant as r —>■ 0. Therefore, as one can see from 
(27) and (28), since r 4 > 0 the sign of f{r) in the vicinity 
of r = 0 depends on the factor 2p—l. Let us discuss the 
cases p > 1/2 and p < 1/2, separately. For p > 1/2, /(r) 
is positive in the vicinity of r = 0. Therefore, one can 
have an extreme black hole with horizon radius r+ = Text 
provided m = rn-ext or q = qext, where 


A4 = — 


(z -|- 1 ) (z -|- 2 ) 
212 ■ 


So, the final form of the metric function is 


(26) 


2p 

^ext 


(z + 2 > 


^r4+2:+2 


r4 


T 4 ,Z ’ 


(32) 


fir) = 1 + 


kP 


rj2p 


pZ -\-2 


r^4-\-z-\-2 ' 


(27) 


^ext — 


(r4 


■2)re,5 , iT4 + z)kPrl^t 


r4 


r4z2 


(33) 


where /c = 0,1 and 


2p ^ (2p - 1) b2(^ 1) 

2Y^l-2p{z-l)-2 




(28) 


Using (16), (22) and (24), it is easy to show that the 
gauge field with the free charge parameter is now given 
by 


_ gib2U-i) 
^r4+i ’ 


It is worthwhile to note that mext and gext are calculated 
from equations fir ext) = 0 and ffext) = 0 , respectively. 
We have also solutions with two inner and outer horizons 
r_ < Text < r^ provided m > mext or q < qext and 
naked singularities when m < mext or q > qext (See Figs. 
(1) and (2)). For p < 1/2, f{r) —>• —00 as r —>• 0 and 
therefore we have Schwarzschild-like solutions. Note that 
in this case £$ < 0 and therefore this result matches with 
the general results we presented before in this section. 
For this case, Figs. (3) and (4) depicts the behavior of 
/(r) with suitable choices of z (see Eq.(31)). 
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FIG. 1: The function /(r) versus r for p = 0.8, z = 3, I — b = 
1, fc = 0 and Text = 1. 


FIG. 4: The function /(r) versus r for p = 1/3 (dashed) 
p = 1/4 (dashdot) p = 0 (solid), z = 9, m = q = l = b= l 
and k = 1. 



FIG. 2: The function /(r) versus r for p = 0.8, z = 3, I — b = 
1, fc = 1 and Text = 1. 



FIG. 3: The function /(r) versus r for p = 1/3 (dashed) 
p = 1/4 (dashdot) p = 0 (solid), z = 9, m = q = l = b= l 
and fc = 0. 


III. HIGHER-DIMENSIONAL LIFSHITZ BLACK 
HOLE/BRANE SOLUTIONS 


In this section, we extend our solutions to the case 
of higher-dimensional solutions. The action of {n -|- 1)- 
dimensional theory (n > 3) can be written as 

^ = -T^ f d"+la;v^|7^-^(V$)2 

loTT I n-l 

-2A -h 


The variation of action (34) leads to the following field 
equations 


|2A -h (2p - 1) 

- V I + 

n — 1 

z=2 ) 

+2^ (iJ,)/ , (34) 

2 = 2 


V^(j) — p'jP 

2=2 

V;. {Hif’') = 0. 


(35) 

(36) 

(37) 


The line element of the higher-dimensional asymptotic 
Lifshitz spacetime is 


ds^ 


12. ^ r^fir) 


+ r^dnl_j_, 


(38) 
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where is an (n— l)-dimensional hypersurface with 

constant curvature (n — l)(n — 2 )k and volume 
Following the method of the previous section, one can 
find the solutions of the field equations (34)-(37) as 


f{r) 

4>(r) 

(^i)* 

(^2), 

(^3)* 


1 - 


m 


kP [n — 2 ) 


o2p 


^n-i+z (2 -f n - 3)^ r2 rr+^+"-i ’ 

(39) 


(n- l)Vz- 
2 

gi62(--i) 




Frr 


z-l 


„n+z — l 


y 2 {n + z — 1) ’ 

^Jk (n — 1) (n — 2) (z — l)r^+"’ 
\/2(z -I- n — 3 )^/ 2 /z- 15 ra -2 


-3 


(40) 

(41) 

(42) 

(43) 


where 


r = z- 2 -f (n-l)/( 2 p-l), 

„2p ^ (2p-l)6^(--i) p 

® (n-l)l- 2 p(^-i)- 2 r’ 


and A = —{z + n — l)(z + n — 2,)/2l'^. Here we pause to 
point out some remarks as in the previous sections. First, 
(^ 3)4 is imaginary in the case of fc = — 1 and therefore 
our solutions rule out the case with hyperbolic horizon. 
Second, At should be finite at infinity since the total mass 
linearly depends on it as we will see later in section (IV). 
Therefore F should be positive. This constraint makes 
some limits on z and p. The allowed ranges for p and 2 
are: 


for p< 1 / 2 , 2 : — 1 > (n — 2 p)/(l — 2 p), 


for 1/2 < p < n/2, all z(> 1) values are allowed, (44) 

for p > n/ 2 , 2 ; — 1 > ( 2 p — n)/( 2 p — 1 ). 

For these allowed ranges ofp and z, the term proportional 
to charge in the metric function disappears at infinity 
and it dominates as t ^ 0. Therefore, since this term 
is positive for p > 1 / 2 , in this case one may encounter 
a black hole with inner and outer horizons or extreme 
black hole for suitable choices of m and q where 


2p 

^ext 


{z + n - 1 ) 

p ' ext 

(n — 2 ) kP 
r(z + n-3) 

(r + ^ + » — 1) z+ra-l 

-p ' ext 


(45) 


(F + z + n — 3) {n — 2)^ kP 
T (z + n — 3)^ 


z+n-3 
' ext 


• ( 46 ) 


On the other hand, for p < 1/2 where charge term is 
negative in metric function, we have Schwarzschild-like 
solutions. 


IV. THERMODYNAMICS OF LIFSHITZ 
BLACK HOLES/BRANES 

In this section, we seek for satisfaction of the first law 
of thermodynamics for our Lifshitz black hole/brane solu¬ 
tions. We start with the calculation of mass that is fun¬ 
damental for thermodynamics discussions. We use the 
modified subtraction method of Brown and York (BY) 
[21]. In order to use the modified BY method [22], the 
metric should be written in the form 

d7?2 

ds^ = -X{n)dt^ + —— + n^dn^. (47) 

I (/vj 

It is obvious from (38) that in our case TZ = r and there¬ 
fore 


Y(7^)=Y(7^) = /(r(7^)). (48) 

The background metric is chosen to be the metric (47) 
with 


Yo(7^) = Yo(7^) = /o(r(7^)) = 1. (49) 


The quasilocal conserved mass can be calculated via 


M = -!- [ d^ipy/a{{Kab-Khab) 

Stt Jb 

(50) 

where a is the determinant of the metric of the boundary 
B, is the extrinsic curvature of the background met¬ 
ric and n“ and are the timelike unit normal vector to 
the boundary B and a timelike Killing vector field on the 
boundary surface, respectively. Using the above modified 
BY formalism, one can calculate the mass of the space 
time per unit volume a;„_i as 


M = 


(n — 1 )to 
167rU+i ’ 


(51) 


where the mass parameter m can be written in term of 
the horizon radius r+ by using the fact that f{r+) = 0 : 


i{r+) = 


, fc/2(n-2)V;+-3 , g2p 


V- ( 52 ) 


(z-l-n —3) ' + 

One can also calculate the charge by using the Gauss law 




'dn. 


( 53 ) 


^ext 
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where and u‘^ are the unit spacelike and timelike nor¬ 
mals to the hypersurface of radius r given as 


= dt = 

y—gtt 


= dr = 
Y 9rr 




Using (53), we obtain the charge per unit volume ujn-i 
as 


Q 


47r 


(54) 


The electric potential U, measured at infinity with re¬ 
spect to horizon is defined by 


U = , (55) 


where x = Cdt is the null generator of the horizon. The 
value of constant C will be explained later. Using (29), 
we can obtain electric potential 


fv^ 


(56) 


The entropy of the black holes can be calculated by using 
the area law of the entropy which is applied to almost all 
kinds of black holes in Einstein gravity including dilaton 
black holes [23, 24]. Thus, the entropy of our solutions 
per unit volume is 


S = 


(57) 



FIG. 5: The behavior of lOT (solid), 10(d^M/dS^)Q (dashed) 
and H^q (dashdot) versus 2 for k = 0 with I = b = 1, 
gi = 0.1, r+ = 0.6, n = 4 and p = 1/3. 


T and electric potential (7 as their conjugate intensive 
quantities, respectively: 



and 



(60) 


Calculations show that intensive quantities calculated by 
(60) coincide with those computed by (56) and (58) pro¬ 
vided C = p. Therefore, these thermodynamics quanti¬ 
ties satisfy the first law of thermodynamics: 


The Hawking temperature can also be obtained as 


dM = TdS + UdQ. (61) 


T. = 




r{r+) 


47rU+i 

1 f (n — 1 -I- z)r^ 


47r 


iz+l 


k{n- 2 f 

(z + n — 3) 


According to (61), it is clear that total mass linearly de¬ 
pends on U. On the other hand U is proportional to At 
^ee Eqs. (55) and (56)). Therefore At should be finite at 
iz+i^^+n-i infinity in order to have a finite mass as we have applied 

Hris condition on our solutions (see Eq. (44)). 

(58) H V 




In order to check the first law of thermodynamics, we 
should write a Smarr-type formula. Using (51), (52), 
(54) and (57), we can write the Smarr-type formula as 


M{S,Q) 


(n-l)(4g)^"~^+^^/^"~^^ 

167rU+^ 

k{n - 1) (n - 2 f 

167rZ^“^ [z + n — 3)^ 

(2p- 1) (45)-^/^"-^) 

7rU-i2(3p-4)/(2p-i)r52(i-2) 

(59) 


We can now consider S and Q as a complete set of exten¬ 
sive quantities for mass M{S,Q) and define temperature 


V. THERMAL STABILITY IN THE 
CANONICAL AND GRAND-CANONICAL 
ENSEMBLES 


In order to study thermal stability of the asymptotic 
Lifshitz dilaton black solutions, we adopt two different 
ensembles, namely, canonical and grand-canonical en¬ 
sembles. In the canonical ensemble, the charge is a fixed 
parameter and therefore the positivity of the heat ca¬ 
pacity Cv = T/{d^M/dS^)Q is sufficient to ensure the 
local stability [25, 26]. This implies that the positivity 
of {d^M/dS‘^)Q guarantees the local stability of the so¬ 
lutions in the ranges where T is positive. It is a matter 






















0.4 


0.2 

0 


- 0.2 


-0.4 



FIG. 6: The behavior of T (solid), 10"^(9^M/aS^)Q (dashed) 
and 10 “^Hs(q (dashdot) versus r+ for A: = 0 with I = b — 1, 
qi = 0.1, 2 = 12, n = 4 and p = 1/3. 



FIG. 7: The behavior of T (solid) and 10 “^H^q (dashdot) 
versus 2 for k = 0 with I — b = 1, qi = 0.2, r+ = 0.55, n = 5 
and p = 2. 


of calculations to show that 

/d^M\ _ z{n-l + 2)r;-"+^ 

\ dS'^ ) Q (n — l)7rP+i 

k (n — 2)^ (z — 

7 rl^~^(n — 1) (z + n — 3) 

2 P(r + n-l)( 2 p- 1) 

+ 7r/(i-2p)(^-i)(n- l)2&2(i-^) ■ 

(62) 

In grand-canonical ensemble Q is no longer fixed. In our 
case, the mass is a function of entropy and charge and 
therefore the system is locally stable provided H^g = 
\d'^M/dSdQ]^ > 0 where the determinant of Hessian ma¬ 


FIG. 8: The behavior of T (solid) and lO'^Hg^g (dashdot) 
versus r+ for k = 0 with t = b — 1, qi = 0.15, 2 = 4, n = 6 
and p = 3. 



FIG. 9: The behavior of lOT (solid), {d'^M/dS^)Q (dashed) 
and (dashdot) versus 2 for k = 1 with I — b = 1, 

qi — 0.05, r+ = 0.6, n = 4 and p = 1/3. 

trix can be calculated as 

JJM ^ p2(n-l + 2)^2p(i-^)-2^;r+.-n+i 

2P-3(2p - l)r62(i--)(n - 
pk (n - 2 f (z - 2)62(^-i)z2p(i-^)r;r-hz-n-i 
2P~^{2p — l)(n — 1) (2 -I- u — 3) 

8p(2p-l)(2-2)<j2^;2“+2 

(n-l)2r&4(i--) ■ 

Let us now discuss thermal stability of the solutions in 
both ensembles. As it is obvious from (62) and (63) 
for 2 = 2 where the allowed range of p is p > 1/2 
(see Eq.(44)), we have thermally stable solutions in both 
canonical and grand-canonical ensembles. For other 
ranges of the parameters, we consider the case k = 0 
and k = 1 separately: 
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FIG. 10: The behavior of T (solid), W~^{d^M/dS^)Q 
(dashed) and (dashdot) versus r+ for fc = 1 with 

I = b = 1, qi = 0.05, z = 12, n = 4 and p = 1/3. 



FIG. 11: The behavior of T (solid), [d'^M / dS^)Q (dashed) 
and 10 “^H;^q (dashdot) versus 2 for k = 1 with I = h — 1, 
qi = 0.1, r+ = 0.6, n = 4 and p = 2. 


i) k = 0: As one can see from Fig. 5, for p < 1/2 the 
behaviors of (^d'^M/dS^) ^ and H^g in terms of 2 are 
opposite. That is in canonical ensemble, there is a Zmax 
such that the solutions are stable provided z < Zmax, 
whereas in the grand-canonical ensemble there is a Zmin 
in which for values greater than it we have thermally 
stable solutions. Also, in this range there is a r+niin in 
the canonical ensemble that for values greater than it 
we have stable solutions while in the grand-canonical en¬ 
semble there is a r_|_ max that for radius of horizon greater 
than it one encounters unstable solutions. This behavior 
is shown in Fig. 6. Of course, one should note that for 
p < 0 where z is always greater than 2 (Eq. (44)), there 
are stable solutions in grand-canonical ensemble. We also 
encounter stable solutions for p > 1/2 for all z(> 1) val¬ 
ues in canonical ensemble while the solutions are stable 



FIG. 12: The behavior of T (solid), W~^{d^M/dS^)Q 
(dashed) and 10“'*H5 (q (dashdot) versus r+ for k = 1 with 
I = b = 1, qi = 0.1, 2 = 1.5, n = 4 and p — 2. 

in the grand-canonical ensemble just for 1 < z < 2. In 
the grand-canonical ensemble when p > 1/2 and z > 2 
there are unstable solutions for z > Zmax > 2 (Fig. 7). In 
latter range there is r_|_ min that for values greater than it 
we have stable black branes (Fig. 8). 

a) k = 1: In this case, solutions are stable in grand- 
canonical ensemble for p < 0 whereas for 0 < p < 1/2 
there are unstable solutions provided z < Zmin (Fig. 9). 
As one can see in Fig. 9, in canonical ensemble, there are 
stable black holes provided z < Zmax when p < 1/2. In 
latter range there is an r_|_ max in grand-canonical ensem¬ 
ble that for r_|_ < r+max black holes are stable while in 
canonical ensemble black holes with horizon radius lower 
than r+min are unstable (Fig. 10). Black holes are al¬ 
ways stable in canonical ensemble if p > 1/2 and z > 2. 
For p > 1/2, we also have stable solutions for 1 < z < 2 
in canonical ensemble provided z > Zmin as it is clear 
from Fig. 11. In grand-canonical ensemble there is also 
a Zmin that for values lower than it solutions are unstable 
for p > 1/2 (Fig. 11). In later range black holes in both 
canonical and grand-canonical ensembles are stable pro¬ 
vided that their radius of horizon r+ > r+ min (See Fig. 
12). 


VI. DYNAMICAL STABILITY OF 
4-DIMENSIONAL ADS BLACK HOLES 

In previous section we investigated the termal stabil¬ 
ity of topological black hole solutions. In addition to 
thermal stability, it is worth discussing stability of black 
hole solutions under dynamical perturbations. There¬ 
fore, in this section, we intend to discuss dynamical sta¬ 
bility of 4-dimensional AdS black hole solutions. It has 
been shown that the properties of perturbations of a 4- 
dimensional static and spherically symmetric background 
under a two-dimensional rotation transformation can be 
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decomposed into odd- and even-parity sectors [27] . These 
perturbations can also be written in terms of the sum of 
spherical harmonics and the dynamical stability is ex¬ 
amined by studying the behavior of perturbation modes. 

Dynamical stability of black hole solutions in Einstein 
gravity coupled to nonlinear electrodynamics (NED) 
sources is presented in [28]. The Hamiltonian of a NED 
Lagrangian C{F) where F = F/A, can be written as 
TLiF) = 2CpF — C. It is convenient to study the dynam¬ 
ical stability in so-called P frame where P = C^pF and 

TLiF) —>• T-L{P). In this frame, there are dynamically sta¬ 
ble solutions under odd-type perturbations provided %p 
has no zero outside the horizon. For even-type perturba¬ 
tions, we have dynamically unstable black hole solutions 
provided 'Hxx > 0 where x = \/—2Q'^ P. In the case of 
AdS solutions (z = 1), one can calculate 


%p 


1 

P 


_4p\ 

) 


(64) 


where P = —^/4. Since 'Hp has no root 
outside the horizon, we have stable solutions under odd- 
type perturbations. Hxx can also be calculated as 


I / 2x‘^ \ 

= -p{2p-l)Q- 

Therefore, we encounter unstable solutions under even- 
type perturbations provided 0 < p < 1/2. 


VII. CLOSING REMARKS 

In this paper, we constructed a new class of asymptoti¬ 
cally Lifshitz topological black hole solutions of Einstein- 
dilaton gravity coupled to a power-law Maxwell field. 
The motivation for taking this kind of Lagrangian orig¬ 
inates from the fact that Maxwell Lagrangian is con¬ 
formally invariant only in four dimensions. In order to 
build a version of Maxwell theory which enjoy the confor¬ 
mally invariance property in all higher dimensions, the 
so-called power-law Maxwell electrodynamics was pro¬ 
posed in Ref. [13]. In fact, in (n -|- l)-dimensional space- 
time the Lagrangian [—is invari¬ 
ant under the conformal transformation 
Afi Afi and 4) —>• 4>, provided p = (n -I- l)/4. We 
first considered the four dimensional action and by using 
this toy model we explained the procedure of obtaining 
the asymptotic Lifshitz topological charged black holes 
of the theory in detail. We allowed the horizon to be a 
hypersurface with k = 1,0,—1 constant curvature. Al¬ 
though we found that the reality of the charge param¬ 
eter rules out the case k = —1 in the presence of dila- 
ton and two electromagnetic fields, a topological Lifshitz 
black hole with fc = — 1 has been obtained in Ref. [4] 


in the presence of a massive electromagnetic field. In¬ 
deed, in the solution obtained in Ref. [4], the massive 
electromagnetic field makes the asymptotic behavior to 
be Lifshitz, while here the set of dilaton and two linear 
electromagnetic fields guarantees the asymptotic Lifshitz 
behavior for topological black hole. We also discussed 
the possible ranges of the parameters for which the met¬ 
ric function goes to /(r) —>• 1 -|- kP/{z^r'^) as r —>• oo 
where z is a dynamical critical exponent shows the de¬ 
grees of anisotropy between space and time in Lifshitz 
spacetime. Then, we generalized our study to all higher 
dimensions. We found that our {n + I)-dimensional solu¬ 
tions are asymptotically Lifshitz provided for p < A12 
and p > n/2, z — 1 > {2p — n)/(2p — 1) while for 
1/2 < p < n/2, all z(> 1) values are allowed. Using the 
modihed Brown York method, we computed the mass of 
the solutions. Next, by computing the thermodynamical 
quantities associated with (n -|- l)-dimensional asymp¬ 
totic Lifshitz dilaton black holes/branes with power-law 
Maxwell field and obtaining the Smarr formula M{S,Q), 
we checked the validity of the first law of thermodynam¬ 
ics, dM = TdS + UdQ, on the horizon. We also explore 
thermal stability of the obtained solutions both in canon¬ 
ical and grand-canonical ensembles. We found that our 
solutions are always stable in both canonical and grand- 
canonical ensembles for z = 2 where the allowed range 
of p is p > 1/2. In grand-canonical ensemble, our black 
hole/brane solutions are always stable for p < 0. In this 
ensemble, there are stable solutions for 0 < p < 1/2 pro¬ 
vided z > Zniin or r < r+niax- For p > 1/2, black branes 
are stable for 1 < z < 2 while they are stable for z < z^ax 
or r > r_|_„iin when z > 2. In latter range of p, black 
holes are stable provided z(r+) > Zniin(?’-i-min)- In canon¬ 
ical ensemble, our black hole/brane solutions encounter 
instability for p < 1/2 and z > Zmax or r < r+min. For 
p > 1/2, black branes {k = 0) are always stable in canon¬ 
ical ensemble while black holes (fc = 1) are stable for 
z > 2. In latter range of p, black hole solutions en¬ 
counter instability in canonical ensemble in the range of 
1 < z < 2 provided z < z„iin or r+ < r_|_„iin- Finally, we 
discussed the stability of 4-dimensional solutions under 
dynamical perturbations in the case of asymptotic AdS 
spacetime where z = 1. 


In this work, we focused on static asymptotic Lifshitz 
topological black holes in the context of Einstein dilaton 
gravity in the presence of electromagnetic fields. This 
work can be extended from different aspects. First, one 
may be interested in studying the effects of rotation on 
thermodynamics and thermal stability of the solutions. 
Second, it is interesting to investigate the thermodynam¬ 
ics of nonlinearly charged topological Lifshitz black holes 
with hyperscaling violation. Also, one may generalize 
these kind of solutions to other kind of nonlinear elec¬ 
tromagnetic fields. Furthermore, this work can be done 
in higher curvature gravities. Some of these works are 
under investigation. 
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